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We report on a search for electromagnetic and/or hadronic showers (cascades) induced by a diffuse 
flux of neutrinos with energies between 5 TeV and 300 TeV from extraterrestrial sources. Cascades 
may be produced by matter interactions of all flavors of neutrinos, and contained cascades have 
better energy resolution and afford better background rejection than through-going ^-induced 
muons. Data taken in 1997 with the AMANDA detector were searched for events with a high- 
energy cascade-like signature. The observed events are consistent with expected backgrounds from 
atmospheric neutrinos and catastrophic energy losses from atmospheric muons. Effective volumes 
for all flavors of neutrinos, which allow the calculation of limits for any neutrino flux model, are 
presented. The limit on cascades from a diffuse flux of v e + z/ M + v T + 77 e + 77 M + z7 T is E 2 44 < 
9.8 x 10 -6 GeVcm -2 s _1 sr _1 , assuming a neutrino flavor flux ratio of 1:1:1 at the detector. The 
limit on cascades from a diffuse flux o£ v e +V e is E 2 ^ <6.5xl0 -6 GeV cm -2 s _1 sr _1 , independent 
of the assumed neutrino flavor flux ratio. 

PACS numbers: 14.60.Lm, 95.85.Ry, 96.40.Tv, 95.55.Vj 



I. INTRODUCTION 

Neutrinos interact principally via the weak force, pos- 
ing a detection challenge for neutrino telescopes but be- 
stowing a valuable advantage on the field of neutrino 
astronomy: neutrino fluxes from astronomical sources 
are essentially unattenuated even over cosmological dis- 



tances. In contrast, high-energy gamma rays are ab- 
sorbed and/or scattered by intervening matter and pho- 
tons, and high-energy cosmic-rays are deflected by galac- 
tic and intergalactic magnetic fields except at the highest 
energies (> 10 19 eV). 

We present a search for the fully-reconstructed light 
patterns created by electromagnetic or hadronic show- 



2 



ers (cascades) resulting from a diffuse flux of high-energy 
extraterrestrial neutrinos. We use data collected in 1997 
from the Antarctic Muon and Neutrino Detector Array 
(AMANDA) for this purpose. Demonstrating ^-induced 
cascade sensitivity is an important step for neutrino as- 
tronomy because the cascade channel probes all neutrino 
flavors, whereas the muon channel is primarily sensitive 
to charged current and T7 M interactions. This is par- 
ticularly relevant in view of the emerging understand- 
ing of neutrino oscillations jy, g, [|, ||, in which the 
flux of Vu would be reduced by oscillations. (The de- 
tection of high-energy atmospheric muon neutrinos by 
AMANDA has been demonstrated by the full reconstruc- 
tion of Cherenkov light patterns produced by up-going 
muons H |(| 0.) Cascades also boast more accurate 
energy measurement and better separation from back- 
ground, although they suffer from worse angular reso- 
lution and reduced effective volume relative to muons. 
Importantly, it is straightforward to calibrate the cas- 
cade response of neutrino telescopes such as AMANDA 
at lower energies through use of, e.g., in-situ light 
sources. Furthermore, cascades become increasingly eas- 
ier to identify and reconstruct as detector volumes get 
larger, so the techniques presented here have relevance 
for future analyses performed at larger detectors. 

Electron neutrinos can produce cascades with no de- 
tectable track via the charged current (CC) interaction 
and all neutrino flavors can produce cascades via the neu- 
tral current (NC) interaction. Cascade-like events are 
also produced in v T CC interactions when the resulting 
r decays into an electron (roughly 18% branching ratio) 
or into mesons (roughly 64% branching ratio) and the r 
energy is below about 100 TeV, at which energy the r de- 
cay length is less than 5 m, so that the shower produced 
by the neutrino interaction and the shower produced by 
the r decay cannot be spatially resolved by AMANDA. 
The contribution of v T to the cascade channel becomes 
important when flavor oscillations are taken into account 
for extraterrestrial |^, [lOj and for atmospheric jllj v- 
induced cascades. For extraterrestrial sources, current 
knowledge of neutrino oscillations suggests a detected 
neutrino flavor flux ratio of i/ e :i/^:i/ r ::l:l:l following an 
expected flux ratio of 1:2:0 at the source. 

The total light output of an electromagnetic cascade 
is approximately 10 8 photons/TeV in ice. Hadronic cas- 
cades have a light yield about 20% lower jL2). An elec- 
tromagnetic cascade develops in a cylinder of about 10- 
15 cm in radius (Moliere radius) and several meters in 
length (about 8.5 m from the vertex of a 100 TeV cascade, 
essentially all charged particles are below the critical en- 
ergy) . Hadronic cascades have longer longitudinal devel- 
opments and larger Moliere radii. As a sparsely instru- 
mented detector, AMANDA is insensitive to the topo- 
logical differences between electromagnetic and hadronic 
cascades. Since the NC interaction has a lower cross sec- 
tion and results in a deposition of less energy than the 
CC interaction, and since we assume a steeply falling 
neutrino energy spectrum, at any given energy a very 



small fraction of the v e events are due to NC interac- 
tions. Hence their impact on the cascade energy resolu- 
tion is small, and the energy spectrum of reconstructed 
cascades closely follows that of the CC v e energy spec- 
trum. 

In this paper, we present limits on the diffuse fluxes of 
(y e + + v T + V e + + v T ) and (y e + V e ), assuming a 
customary E~ 2 power law spectrum at the source. These 
limits are based on the observation of no events consistent 
with a diffuse flux of high-energy extraterrestrial neutri- 
nos. We also present effective volumes for all neutrino 
flavors to facilitate the calculation of a limit for any flux 
model. (A search for up-going muons produced by a 
extraterrestrial diffuse flux is presently being conducted 
and preliminary results have been reported in JOI.) 



II. THE AMANDA-BIO DETECTOR 

The data used in this work was taken with the 
AMANDA-B10 detector in 1997. AMANDA-B10 §, @, 
|l4| was commissioned in 1997 with a total of 302 op- 
tical modules (OMs) arranged on 10 strings, at depths 
between 1500 m and 2000 m below the surface of the ice 
at the South Pole. The strings are arranged in two con- 
centric circles 35 m and 60 m in radius, with one string at 
the center. The OMs in the inner four (outer six) strings 
have a 20 m (10 m) vertical separation. Each OM con- 
tains a 20 cm photo-multiplier tube (PMT) in a spherical 
pressure vessel. Coaxial cables in the inner four strings 
and twisted pair cables in the outer six strings provide 
high voltage to the PMTs and simultaneously transmit 
their signals to the electronics housed on the surface. The 
detector is triggered using a majority condition in which 
an event is recorded if more than 16 modules have a sig- 
nal (i.e., were "hit") in a 2 /us time window. A total of 
1.05 x 10 9 events were recorded during an effective live- 
time of 130.1 days. 

The optical properties of the ice have been studied with 
in-situ light sources and with atmospheric muons. These 
studies have shown that ice at the South Pole is not per- 
fectly homogeneous, but rather consists of horizontal lay- 
ers corresponding to global climatological conditions in 
the past, such as ice ages. These layers lead to a modu- 
lation of the absorption and effective scattering lengths 
as a function of depth Jl[|. Optical properties are also 
modified by the presence of drill-hole bubbles which are 
created during the drilling and deployment processes. 



III. METHODS FOR CASCADE 
RECONSTRUCTION 

Simple reconstruction algorithms are initially applied 
to the data. These methods are used to reduce the data 
sample size and to seed more sophisticated reconstruc- 
tion algorithms, while maintaining high passing rates for 
simulated signal events. For cascades, the mean posi- 
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tion of the hit OMs, or center of gravity, is used as the 
first guess of the position. In order to efficiently reject 
muons, they too are reconstructed, beginning; with a first 
guess track fit called the line fit |L(|. The line fit is an 
algorithm that assumes that hits can be projected onto a 
line, and that the particle which produced the hits trav- 
els with a velocity vn ne and has a starting point r$. The 
fit minimizes the quantity ^^ ta (fj — f — vn nc ■ ti) 2 as 
a function of r*o and vn no , where Nhits is the number of 
hits in the event. These procedures are described in more 
detail elsewhere @ Q . 

After calculating the first guesses, three maximum like- 
lihood methods are used consecutively to reconstruct pre- 
cisely the cascade vertex position, time, energy and di- 
rection. These methods are described below. 



A. Single Photoelectron Vertex Position and Time 
Reconstruction 

The cascade vertex position and creation time are re- 
constructed using a maximum likelihood function that 
takes into account the Cherenkov emission, absorption 
and scattering of light. This vertex information is re- 
quired for rejecting potential backgrounds and for sub- 
sequent fits for energy and direction. This procedure is 
quite similar to the algorithms used for muon fitting [jl7| . 
A more comprehensive description of the different cas- 
cade reconstruction methods can be found in . 

We use a likelihood function having the form: 

N hits 

u.p(tL,di), (i) 

i=0 

where t rcs = t^it — ^Cher is the difference between ob- 
served hit time and expected time for Cherenkov emis- 
sion without scattering-the time residual-and p(t Tes , d) 
is the probability of observing a photon at a time resid- 
ual i res at a distance d from the emitter. The label "spe" 
indicates that £l p t c assumes all hits are due to single pho- 
toelectrons. 

The probability p(t TCS , d) was generated by parametriz- 
ing simulations of light propagation in ice. The product 
in Eq. |l| is calculated using all hit OMs. The maximiza- 
tion of £t p t e provides a good estimate of the vertex posi- 
tion and time of the cascade. 



B. Multi-photoelectron Vertex Position and Time 
Reconstruction 

The single photoelectron likelihood can be refined by 
taking into account that the time measured in each PMT 
is the time of the first photon to be observed. If a PMT 
receives N photons, the probability of measuring a time 
residual, i res , and the associated likelihood function are: 



(roc \ N — l 

jf dt'p(t',d)J (2) 

JV hit3 

£™r= Hp(N,tie S ,di) (3) 

where the "mpe" label indicates that £~ lpe describes 
multi-photoelectron hits. In AMANDA we use the mea- 
sured pulse amplitude as an estimator for N. 

The maximization of £™ pe is used to estimate the most 
likely vertex position x and time t of a cascade. The 
multi-photoelectron vertex reconstruction uses the max- 
imization of C S g° to seed initial values of cascade vertex 
position and time. 

C. Energy and Direction Reconstruction 

Cascade energy and direction are reconstructed using 
a likelihood function assembled from the probabilities of 
an OM being hit or remaining unhit assuming a cascade 
hypothesis: 

Hit OMs 

£e,A = JJ P Hit(rfi,^,^) X 
i=0 
Unhit OMs 

Pvnbit{di,E,n) (4) 

i=0 

Maximization of Cea provides the most likely value of 
energy E and direction n of a cascade. Note that in 
principle this procedure also allows for the reconstruc- 
tion of position (but not time) of a cascade. Monte Carlo 
studies have shown, however, that the position resolution 
obtained by maximizing Lea is not as good as that ob- 
tained with £l p 4 ° or £™ pc . 

IV. DETECTOR RESPONSE TO CASCADES 

In the absence of a tagged source of high-energy 
neutrino-induced cascades, to understand the response 
of the detector we rely on in- situ light sources, catas- 
trophic energy losses by down-going cosmic-ray muons, 
and Monte Carlo simulations. The successful reconstruc- 
tion of these data demonstrate detector sensitivity to cas- 
cade signals. 

A. Pulsed Laser 

A pulsed laser operating at 532 nm on the surface is 
used to send light through optical fibers to diffuser balls 
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FIG. 1: The plots show results of the reconstruction of pulsed 
laser data (solid) and simulation (dashed) with the diffuscr 
ball next to OM 69 (situated near the center of the detector). 
The differences between simulated and reconstructed x and 
z components of the position are shown on top. The speed 
of the line fit, wii Ile , is shown in the bottom left. The bottom 
right plot shows the reduced L™^ c , described in the text. The 
vertical scale for all four plots is arbitrary. The position of the 
OMs is known with about 1 m precision [ pl| , so there is no 
discrepancy in the plot in the upper right. The discrepancy 
in the plot in the lower right arises due to the simplified ice 
model used for the pulsed Laser. 



embedded in the ice close to almost every OM in the 
detector. A comparison of Monte Carlo and experimen- 
tal data for these in-situ light sources deepens our un- 
derstanding of reconstruction performance and detector 
signal sensitivity. The photon intensity that can be pro- 
duced at each diffuser ball is not known a priori, so we 
force the number of hit channels in experimental and sim- 
ulated pulsed laser data to match. Thus, the simulations 
predict that the laser produces pulses in the ice com- 
prising 5 x 10 7 — 1 x 10 9 photons (corresponding to a 
maximum cascade energy of roughly 10 TeV). The laser 
pulses are roughly 10 ns wide, short enough to mimic the 
time structure of true cascades. 

Although highly useful as a cascade calibration source, 
the pulsed laser system has some minor drawbacks. The 
diffuser ball light output is expected to be isotropic, so 
the laser data does not provide information about the an- 
gular response of the detector to cascades. The laser pro- 
duces light at A = 532 nm and at this wavelength the op- 
tical ice properties are different from those at Cherenkov 
radiation wavelengths. The effective scattering length at 
532 nm is 18-30 m and depends on depth. The absorption 
length at 532 nm is 25 m and independent of depth 
(at the shorter wavelengths characteristic of Cherenkov 



radiation the absorption length is about 100 m). 

Independent data sets taken with diffuser balls in a va- 
riety of locations are reconstructed with the first guesses 
and with the time-position reconstruction algorithm de- 
scribed above. The position resolution is about 1 m in 
the z dimension and about 2 m the in x and y. It is 
better in z due to closer OM spacing in that dimension. 

The pulsed laser simulation uses a simplified optical 
model of ice properties: the drill-hole bubbles are taken 
into account, but no depth dependence is used for the 
scattering length. In spite of this simplification, the ver- 
tex resolution of the pulsed laser data agrees well with 
simulations, showing that the detector can be used to 
reconstruct the position of contained point-like events. 
(Contained events are defined as events whose recon- 
structed vertex lies within a right cylinder of height 400 m 
and radius 60 m, centered on the AMANDA-B10 detec- 
tor and encapsulated by it.) Figure [l] shows the results 
of the reconstruction of pulsed laser data and simulation. 



B. Catastrophic Muon Energy Losses 

The vast majority of the events recorded by AMANDA 
are down-going muons induced by cosmic-ray air showers. 
This background has been simulated with the Corsika 
program p2| using the average winter air density profile 
at the South Pole and the QGSJET hadronic model |25( | 
option. The cosmic-ray composition is taken from p4[ . 
The propagation of muons through ice is simulated with 
the program Mudedx [^5[ ^6). Optical properties of the 
ice, including depth dependence and drill-hole bubbles, 
are also simulated. 

From the large sample of atmospheric muons it is 
possible to extract a subset of cascade-like events in 
which the majority of the recorded hits come from 
catastrophic, localized energy loss of the muon (e.g., a 
bright bremsstrahlung) . The extraction of these events 
is achieved using criteria 1-9 from table ||, i.e., we do not 
require these events to reconstruct as up-going cascades. 
(We do, however, still reject obvious down-going muons 
via criterion number 8.) Based on the number of hits 
produced by the brightest cascade in simulated events, 
and on a visual study of these events, we have confirmed 
that after applying these selection criteria the remaining 
events are indeed cascade-like. 

Figure || shows the energy spectra of muon energy 
losses for experimental data and simulated atmospheric 
muons. The experimental and simulated data agree rea- 
sonably well, but not perfectly. The difference is dis- 
cussed in Section VI. 



C. Monte Carlo Prediction for Neutrino-Induced 
Cascade Reconstruction 

To study the performance of the reconstruction al- 
gorithms we simulated a flux of v e + T> e following an 
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FIG. 2: Energy spectra of reconstructed atmospheric muon 
energy losses for experimental data and standard simulated 
muon background after application of selection criteria 1-9 
from table |. Agreement between simulation and experimen- 
tal data depends on a combination of the simulated ice prop- 
erties and the OM angular and absolute sensitivity. These 
effects have been taken into account in the calculation of sys- 
tematic uncertainties in Section |vi|. 

E~ 2 power law spectrum. Neutrinos from astrophysical 
sources are expected to have a hard spectrum, reflect- 
ing the processes in the cosmic accelerators that generate 
them. Earth absorption and NC scattering are taken into 
account in the simulation. The "Preliminary Reference 
Earth Model" is used to calculate the Earth's density 
profile . We calculate differential cross sections using 
CTEQ5 following Gandhi et al. p8| . For v T interactions, 
the simulation of the r decay uses TAU0LA [3(| Q. 

Neutrino-induced cascades are reconstructed following 
the procedure described in Section III. Position, zenith 
angle and energy resolutions for a flux of v e + 77 e are cal- 
culated using the difference distributions shown in Fig. [|. 
The position resolution is roughly 4 m in the z dimen- 
sion and 5.5 mini and y for contained cascades. (Note 
that position resolution obtained with the pulsed laser is 
better than that predicted for Cherenkov light because 
optical ice properties are more favorable at the longer 
wavelength.) The reconstructed position is biased in the 
direction of the cascade, but since the mean of this shift 
is only about 2 m for contained cascades it has a negligi- 
ble impact on the final result. Zenith angle resolution is 
25°-30° depending on the cascade energy. 

The energy reconstruction has a resolution of 0.12- 
0.20 in log 10 -E for contained cascades in the range 1 TeV 
- 100 TeV, increasing as a function of cascade energy. En- 
ergy reconstruction of contained cascades is possible from 
approximately 50 GeV (the minimum energy cascade 



FIG. 3: The four plots show the difference between simulated 
and reconstructed vertex position, energy and direction of 
cascades. The monoenergetic cascades have E — 10 TeV and 
are contained within the detector. Vertical scale for all four 
plots is arbitrary. 

which can trigger the detector) to about 100 TeV. At en- 
ergies higher than 100 TeV all, or almost all, of the OMs 
are hit and thus energy reconstruction by the minimiza- 
tion of Eq. ^ is not possible. (Such high-energy events 
would, however, certainly be identifiable, and probably 
they can be reconstructed by other techniques.) 



V. ANALYSIS 
A. Filter 

The first step in the analysis is to apply an initial set 
of selection criteria, here called the "filter," which results 
in a reduction of the data sample size by more than two 
orders of magnitude. The filter first removes spurious hits 
arising from electronic and PMT noise. It then uses the 
fast reconstruction algorithms described earlier, a simple 
energy estimator, and topological characteristics of the 
hit pattern to select potential signal events. The various 
filter steps were tuned to reject a simulated background 
of down-going cosmic-ray muons. 

After the filter has been applied, hits likely to have 
come from cross-talk are removed, as explained below. 
Then each event is reconstructed twice, first with a cas- 
cade hypothesis and then with a muon hypothesis. Sev- 
eral selection criteria are then applied to the data based 
on the results of the reconstruction. 

Atmospheric v e and neutrinos are simulated accord- 
ing to the flux calculated by Lipari p2| . Contributions to 
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correlation in time and amplitude of real and cross-talk 
hits. The map shows that cross-talk occurs for OMs in 
the same string that are close neighbors or for OMs in the 
same string that are separated by 150-200 m. The origin 
of cross-talk is correlated with the relative positioning of 
individual electrical cables within the string |33| . 

Cross-talk hits may also be characterized by narrow- 
ness (small time over threshold or "TOT") coupled with 
unexpectedly large amplitude. Cross-talk-induced and 
light-induced hits lie in different regions of the amplitude 
vs. TOT space and may therefore be separated from one 
another. Figure ^ shows the amplitude vs. TOT distri- 
butions for both types of hits. 

The cross-talk map and the amplitude vs. TOT in- 
formation are both used to remove cross-talk from the 
experimental data. 



Selection Criteria 



FIG. 4: Amplitude vs. time over threshold (TOT) distri- 
bution for hits on OM 149 due to light and cross-talk. The 
data shown to the right of the solid curve are generated using 
the pulsed laser with only the high voltage for the PMT in 
OM 149 enabled (all other PMTs had their high voltage dis- 
abled) . This region of the plot therefore contains hits created 
by light. The data shown to the left of the solid curve are 
also generated using the pulsed laser, but with a diffuser ball 
located 200 m above OM 149, close to OM 129, in the same 
string as OM 149. Only OM 129 had its high voltage enabled. 
This region of the plot therefore contains hits in OM 149 due 
to cross-talk. Hits are removed from an event if they lie to 
the left of the dashed curve or if they have a TOT smaller 
than 125 ns (indicated by the dashed line). 

cascades from both atmospheric v e (CC and NC interac- 
tions) and (NC interactions) are taken into account. 
In the energy range relevant to this analysis (E >5 TeV), 
neutrino oscillations are not important in the simulation 
of atmospheric neutrinos. 



B. Removal of Electronic Cross- Talk 

Electronic cross-talk is present in the twisted pair ca- 
bles used for strings 5-10. Spurious hits arising from 
cross-talk can degrade the reconstruction quality. Cross- 
talk is not included in the simulations, so its removal is 
an important facet of this analysis. 

We generate a detector-wide map of the cross-talk us- 
ing the pulsed laser. For each OM in strings 5-10, pulsed 
laser data are taken in which only the PMT in the OM 
near the laser diffuser ball has its high voltage enabled. 
Any hit in the flashed OM is thus known to be due to light 
and any hit in any other OM is known to be due to cross- 
talk. The cross-talk map identifies the pairs of cross- 
talk-inducing and cross-talk-induced OMs as well as the 



Table | lists the selection criteria and the passing rates 
for experimental data and the various samples of Monte 
Carlo used in this analysis. Selection criteria which have 
not already been described in Section III are described 
below, followed by a physical justification for each crite- 
rion. 

The ratio of the smallest to the largest eigenvalues of 
the tensor of inertia of the position of the hits, A1/A3 
or sphericity, is used to classify events 1 . Small values of 
the sphericity correspond to hits located along a narrow 
cylinder, as expected for a muon. Values of the sphericity 
close to unity correspond to a spherical distribution of the 
hits, as expected for contained cascades. 

The dispersive nature of the ice and of the cables that 
transmit the electrical signals from the OMs to the elec- 
tronics on the surface render sharp signals in the ice 
into significantly broader pulses at the surface. For this 
reason, counting the number of OMs with large TOT, 
-WbigTOT, gives a rough estimate of the energy of the 
event. For this analysis a TOT is considered large if it 
exceeds the value estimated from Monte Carlo simula- 
tions to correspond to one photoelectron by at least a 
factor of 1.5. A contained cascade with 300 GeV of en- 
ergy corresponds roughly to iVbi g TOT = 6. 

The quality of the likelihood reconstruction is deter- 
mined from the reduced likelihood parameter, defined by 
L = —logC/{Nhxt, s — Nftt), where N& t is the number of fit- 
ted parameters. Lower values of L correspond to better 
reconstruction quality. 

A hit is considered direct if the time residual is between 
-15 ns and 75 ns. The number of direct hits, JVdi r , is an- 
other measure of the quality of the reconstruction. Both 



1 To calculate a tensor of inertia in this context, a unit mass is 
hypothesized at each hit OM position. 
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N™ pc > 12 
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1.26 10 6 
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32.1 


26.2 
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2520 


8 


M > 80° 


3.62 10 5 


2.47 10 4 


22.5 


18.5 
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1676 


9 


Slices in z c 


1.48 10 5 


1.19 10 4 


10.6 


8.5 


528 


711 


10 


cos(6» c ) < -0.6 


675 


84 


1.3 


1.0 


106 


156 


11 


E c vs. p m P c 








0.01 


0.01 


28.7 


43 



TABLE I: Selection criteria used in the search for high-energy z/-induced cascades with AMANDA-B10. The number of events left after 
applying each selection criterion to experimental data and background simulations of atmospheric fi, u e and are shown. We simulated 
20.3 days of atmospheric fi data, and 130.1 days of atmospheric v e and v^. Signal simulation is also shown for i>; +77; and v e +F e assuming 
E~ 2 spectra, a flux of 1 X 10~ 4 GeV^ 1 cm~ 2 s _1 sr _1 . 



the single and multi-photoelectron position-time recon- 
struction report the number of direct hits. 

Criteria 1-5 of table [| correspond to the filter and these 
criteria must be satisfied by all events in the analysis. 
The cut on iVbigTOT selects high-energy events, and the 
cut on the sphericity selects events in which the hit topol- 
ogy is cascade-like. The cut on the speed of the line fit, 
viine, removes easily identified down- going muons. As 
shown in Fig. || and Fig. ||, cuts on the likelihood pa- 
rameter (criterion 6) and the number of direct hits (cri- 
terion 7) are used to eliminate non-cascade-like events 
and preserve cascade-like events with good reconstruc- 
tion quality. In addition to the filter criteria, angular 
cuts (criteria 8 and 10) are used to reduce clear muon- 
like events and to select up-going cascade-like events, and 
criterion 11 selects high-energy events within a given dis- 
tance from the vertical axis of the detector, where E c is 
the recon structed cascade en ergy using Eqs. § and |J, and 
p mpc = y/(x m P c ) 2 + (y m P°) 2 . (Criterion 9 is discussed in 
detail below.) 

The filter was developed based strictly on the predic- 
tions of the signal and background Monte Carlo. As more 
and more cuts were applied after the filter, it was found 
that the experimental data and simulations disagreed in 
the shape of the z-component of the reconstructed cas- 
cade position and in the reconstructed cascade direction. 
Inadequately simulated or unsimulated detector instru- 
mental, such as optical properties of the ice and cross- 
talk, contribute to this disagreement. Restricting the re- 
gions of the detector used in this analysis reduces the ef- 
fective volume by more than a factor of two, but restores 
the agreement between experimental data and simula- 
tions. Moreover, the disagreement is also present in the 
reconstructed cascade direction. Only the regions of the 
detector that are accepted, show agreement in the re- 
constructed cascade direction. Events are accepted only 
if their reconstructed vertices satisfy selection criterion 
number 9: —80 m < z c < —40 m or 40 ru < z c < 160 m 



0.3 



0.2 



0.1 







Experiment 

fi background MC 

v Signal MC 

u Signal MC — Contained 



5 5.5 6 6.5 7 7.5 8 8.5 9 
Likelihood Parameter L mpe 



FIG. 5: Normalized distribution of the reduced likelihood pa- 
rameter L?? e for experimental data (solid line), background 
atmospheric fi simulation (dashed line), v e simulation assum- 
ing an E~ 2 power law spectrum (dotted line) and contained 
Ve simulation assuming an E~ 2 power law spectrum (dotted- 
dashed line). Contained events have their vertices in a cylin- 
der 400 m in height and 60 m in radius, roughly matching the 
detector dimensions. Selection criteria 1-5 from the same ta- 
ble have already been applied to all the samples shown. The 
arrow indicates the region removed by cut 6 from table |. 



with respect to the center of the detector (located at 
1730 m below the surface). 
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FIG. 6: Normalized distribution of direct hits for experi- 
mental data (solid line), background atmospheric /i simula- 
tion (dashed line), v e simulation assuming an E~ 2 power law 
spectrum (dotted line) and contained v e simulation assuming 
an E~ 2 power law spectrum (dotted-dashed line). Contained 
events have their vertices in a cylinder 400 m in height and 
60 m in radius, roughly matching the detector dimensions. 
Selection criteria 1-5 from the same table have already been 
applied to all the samples shown. The arrow indicates the 
region removed by cut 7 from table |l| 



VI. SYSTEMATIC UNCERTAINTIES 

There are several uncertainties inherent in estimating 
the detector sensitivity to high-energy neutrino-induced 
cascades. First, the detection medium is a natural ma- 
terial (South Pole ice) whose properties are not pre- 
cisely known. Second, there are no sufficiently powerful 
accelerator-based sources of neutrinos available for use as 
calibration beams. Consequently, the understanding of 
the detector sensitivity is achieved using down-going at- 
mospheric muons, in-situ light sources and Monte Carlo 
simulations. 

To estimate the systematic uncertainty due to impre- 
cise knowledge of the optical properties of the ice, simu- 
lations have been performed using the least and the most 
transparent ice that we have measured at AMANDA 
depths. The cascade sensitivity is modified by 20% us- 
ing either extreme model of the optical properties. Un- 
certainties in the bubble density in the drill-hole ice 
translate to uncertainties in the OM angular sensitivity. 
Monte Carlo simulations with increased bubble density 
in the drill-hole ice degrade the cascade sensitivity by 9%. 
The absolute sensitivity of the OMs is also uncertain at 
the level of 40%. Monte Carlo simulations with altered 
absolute OM sensitivity modifies the cascade sensitivity 



by 5%. This dependence is weak due to the high Nhits 
requirement imposed by this analysis. (The dependence 
is much stronger at earlier stages of the analysis, where 
the average Nhits is much lower. For example, before the 
filter is applied a variation in absolute OM sensitivity of 
40% results in a modification of the cascade sensitivity 
by roughly 35%.) 

Cross-talk can reduce the sensitivity of the detector 
to high-energy neutrino-induced cascades. Events for 
which cross-talk is not fully removed are typically mis- 
reconstructed and are therefore unlikely to have suffi- 
cient quality to pass our selection criteria. The pulsed 
laser data is used to estimate the cascade sensitivity loss 
due to cross-talk for different locations in the detector. 
These studies indicate that the sensitivity is degraded by 
7% due to cross-talk. (This 7% degradation is applied 
directly to the limit and not treated as a systematic un- 
certainty.) Related to cross-talk is the uncertainty in the 
limits due to using slices in z c . Changing each boundary 
of the slices by the position resolution in z modifies the 
cascade sensitivity by 4%. 

Uncertainties in the limits due to neutrino-nucleon 
cross sections, total cascade light output, and cascade 
longitudinal development have also been estimated us- 
ing Monte Carlo simulations. For each of these cases the 
cascade sensitivity is modified by < 5%. 

The systematic uncertainties discussed so far are added 
in quadrature, giving an overall systematic uncertainty 
on the sensitivity of 25%. We follow the procedure de- 
scribed in [55| to determine how to modify the final 
limit in light of this systematic uncertainty, assuming 
that the uncertainties are of a Gaussian nature. 

The spectrum of cascade-like events produced by 
dow n-goin g muons is shown in Fig. || (see also Sec- 



tion IV B ). Standard simulations as well as simulations 
with modified ice properties and OM angular and abso- 
lute sensitivities have been performed. The disagreement 
between experiment and simulations may be explained by 
the uncertainties in the knowledge of the optical proper- 
ties of ice, the OM sensitivity, the cosmic-ray spectrum 
and the rate of muon energy losses. From Fig. |^ it can be 
seen that reasonable agreement between experiment and 
simulations is restored by shifting the energy scale by up 
to 0.2 in log 10 E. This uncertainty in the energy scale 
results in an uncertainty on the sensitivity of less than 
25%. This uncertainty is not independent of the other 
sources of systematic uncertainty that we have studied. 
It demonstrates, however, that the overall systematic un- 
certainty has not been grossly under- or overestimated. 



VII. RESULTS 

The analysis is applied to simulated samples of at- 
mospheric v e and background, high-energy neutrino 
signal (all flavors), and atmospheric muons, and to the 
1997 experimental data set. In the experimental data 
zero events are found. The simulation of atmospheric v e 
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predicts 0.01 events, and the simulation of atmospheric 
i/fj, predicts 0.01 events from NC interactions (both these 
numbers have been rounded up from distinct smaller val- 
ues). Zero events are found in the simulated atmospheric 



muon sample after all cuts. A limit on the flux of neu- 
trinos assuming an E~ 2 power law spectrum is set using 
the following formula: 



J 



E 



dE 



N< 



I 



\dfldE 



(5) 



where I is the neutrino flavor, E the neutrino energy, 
9 the neutrino zenith angle, Ng % = 2.62 determined 
using the unified Feldman-Cousins procedure ]36| with 
a correction applied for the estimated 25% systematic 
uncertainty @ ||, T the live-time (130.1 days), N A 
Avogadro's number, p- lcc the density of ice, a\ ot (E) the 
neutrino cross section p8| , V^q(E, 8) the effective vol- 
ume of the detector (see table ||), // the fraction of the 
total neutrino flux comprised by the neutrino flavor I, 
and £i(E,6) a function that corrects the flux for Earth 
absorption and NC scattering. The integration of Eq. || 
has been done for neutrino energies between 5 TeV and 
300 TeV. 

The 90% C.L. limit on the diffuse flux of v e + + v T + 
v e + 77,, + V T for neutrino energies between 5 TeV and 
300 TeV, assuming a neutrino flux ratio of 1:1:1 at the 
detector, is: 



E 2 ^- < 9.8 x 10~ 6 GeVcm^s^si-- 1 . (6) 
dE 

The 90% C.L. limit on the diffuse flux of v e + 77 e for 
neutrino energies between 5 TeV and 300 TeV is: 



E 2 ^- < 6.5 x 10~ 
dE 



GeVcm 2 s 1 sr 



(7) 



The latter limit is independent of the assumed 
neutrino flux ratio. The limits without in- 
corporating the effects of systematic uncer- 
tainties are 9.1 x 10 -6 GeV cm -2 s _1 sr -1 and 
6.1 x 10~ 6 GeV cm~ 2 s _1 sr -1 , respectively. (Note 
that since the limit in Eq. ^| is on the sum of the fluxes 
of all neutrino flavors, and the limit in Eq. [?] is on an 
individual flavor, the former limit should be divided by 
a factor of three to compare it properly to the latter.) 

Our results together with other limits on the flux of dif- 
fuse neutrinos are shown in Fig. |^. Since recent results 
from other low energy neutrino experiments |jL ^, |j| |^] in- 
dicate that high-energy cosmological neutrinos will have 
a neutrino flavor flux ratio of 1:1:1 upon detection, in this 
figure we scale limits derived under different assumptions 
accordingly. For example, to do a side- by-side compari- 
son of a limit on the flux of v e + -\-v T + V e -\-V fJi +V T , 
derived under the assumption of a ratio of 1:1:1, to a limit 
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FIG. 7: Distribution of and v T energies after all se- 

lection criteria have been applied. The relative normaliza- 
tion between the histograms indicates the relative number of 
events for each neutrino flavor that passes all the selection 
criteria. The initial energy distribution follows an E~ 2 spec- 
trum. Neutrino absorption inside Earth, NC scattering and 
r de cay h ave been taken into account as described in Sec- 
tion live. 



on just the flux of Vr, 



the latter must be degraded 



by a factor of three. (N.B.: We assume that v:V::l:\.) 
Following the Learned and Mannheim prescription for 
presenting limits |37j] , we show neutrino energy distribu- 
tions after applying all the selection criteria in Fig. |7[ 



It should be noted that most searches of diffuse fluxes 
shown in Fig. |^ are based on the observation of up-going 
neutrino-induced muons. Only Baikal and AMANDA 
have presented limits from analyses that search for 
neutrino-induced cascades and only the AMANDA anal- 
ysis uses full cascade event reconstruction. 
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Baikal (v c )[6] 
(prelim.) 
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(prelim.) 

VMANDA-BlOt^+Vj) 

Baikal NT96+NT200 (Vj+v,) (prelim.) 
^AMANDA-BIO (V +V ) (prelim.)[3] 



[X] Multiplicative factor of X 
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not assume 1:1:1 flux ratio 
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FIG. 8: The limits on the cascade-producing neutrino 
flux, summed over the three active flavors, presented in 
this work and in other experiments, with multiplicative fac- 
tors applied as indicated to permit comparison of limits de- 
rived with different assumed neutrino fluxes at the detec- 
tor: Baikal t p e ) |3£| ] (at the resonance); Baikal NT96 
(^M+F^+^e) f§; Frejus K+77 M ) Ggl; MACRO K+F M ) 0. 
Baikal NT96+NT200 (v, + Vi) [§|, AMANDA B-10 

(y^ + Jig) ); Also shown are the predicted horizontal and 
vertical v e and atmospheric fluxes 



VIII. CONCLUSIONS 

High-energy neutrino-induced cascades have been 
searched for in the data collected by AMANDA-BIO in 



1997. Detailed event reconstruction was performed. Us- 
ing in-situ light sources and atmospheric muon catas- 
trophic energy losses, the sensitivity of the detector to 
high-energy cascades has been demonstrated. 

No evidence for the existence of a diffuse flux of neu- 
trinos producing cascade signatures has been found. Ef- 
fective volumes as a function of energy and zenith an- 
gle for all neutrino flavors have been presented. The ef- 
fective volumes allow the calculation of limits for any 
predicted neutrino flux model. The limit on cascades 
from a diffuse flux of v e + v a + v T + v e + 77„ + v T is 



cm 2 s _1 sr _1 assuming a neu- 
trino flavor flux ratio of 1:1:1 at the detector. The limit 



E 2 ^ < 9.8 x 10~ 6 GeVc™ 2 - 1 - 1 



on cascades from a diffuse flux of v e 



is < 



6.5 x 10 GeVcm s sr , independent of the as- 
sumed neutrino flux ratio. The limits are valid for neu- 
trino fluxes in the energy range of 5 TeV to 300 TeV. 
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TABLE II: Effective volume, in units of 10 3 km 3 , for all neutrino flavors as a function of energy and zenith angle after all the 
selection criteria have been applied. Uncertainties are statistical only. 
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